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Abstract 

The planewave response of a linear passive material generally cannot be characterized by a single scalar 
refractive index, as directionality of energy flow and multiple wavevectors may need to be considered. 
This is especially significant for materials which support negative refraction, negative phase velocity, and 
counterposition. By means of a numerical example based on a commonly studied bianisotropic material, 
our theoretical investigation revealed that: (i) negative (positive) refraction can arise even though the 
phase velocity is positive (negative); (ii) counterposition can arise in instances of positive and negative 
refraction; (iii) the phase velocity and time-averaged Poynting vectors can be mutually orthogonal; and 
(iv) whether or not negative refraction occurs can depend upon the state of polarization and angle of 
incidence. A further numerical example revealed that negative phase velocity and positive refraction can 
co-exist even in a simple isotropic dielectric material. 

1 Introduction 

Negatively refracting materials have been subjected to intense interest since the first experimental demon- 
stration of a negatively refracting isotropic dielectric-magnetic material at the turn of this century [1, 2]. 
The enhanced scope for negative refraction in anisotropic materials [3, 4], and in materials which exhibit 
either isotropic [5, 6, 7] or anisotropic [8] magnetoelectric coupling, has been widely reported upon. In- 
deed, bianisotropy may be responsible for unexplained features in certain materials which were taken to be 
isotropic dielectric-magnetic materials [9]. 

We confine ourselves here to linear passive materials. The characterization of negative refraction in 
anisotropic and bianisotropic materials should be expected to be much more complicated than in isotropic 
dielectric-magnetic materials, as the effects of directionality, magnetoelectric coupling, and two refraction 
wavevectors need to be considered [10]. However, it is often overlooked that even in an isotropic dielectric 
material the effects of directionality can be important if nonuniform plane waves are considered. Negative 
phase velocity (NPV) — which means that the phase velocity of a plane wave casts a negative projection 
onto the time-averaged Poynting vector — is commonly taken as convenient indication of the propensity 
for negative refraction [11, 12]. For uniform plane waves in isotropic dielectric-magnetic materials, the 
phase velocity and time-averaged Poynting vector are either parallel or antiparallel. Accordingly, as regards 
uniform planewave propagation, negative refraction and NPV are held to be effectively synonymous terms 
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for these materials. But for nonuniform plane waves, the time-averaged Poynting vector and the phase 
velocity vector are not necessarily parallel or antiparallel, even in isotropic dielectric materials [13]. The 
introduction of anisotropy or bianisotropy further complicates the issue, with the phase velocity and time- 
averaged Poynting vector being generally neither parallel nor antiparallel, for both uniform and nonuniform 
plane waves. Thus, NPV should not be generally assumed to be a definite signature of the capability to 
exhibit negative refraction. 

For many practical applications, the direction of energy flow, as delineated by the time-averaged Poynting 
vector, and its deflection at the planar boundary between two different mediums may be more significant 
than the deflection of the wavevector. It is quite possible for the real part of a refraction wavevector and its 
associated time-averaged Poynting vector to be oriented on opposite sides of the normal to a planar interface. 
This counterposition of the real part of the refraction wavevector and the time-averaged Poynting vector 
has been theoretically demonstrated as taking place in certain anisotropic [14, 15, 16] and bianisotropic [17] 
materials. Furthermore, counterposition can also arise for nonuniform planewave propagation in certain 
isotropic dielectric materials, as we describe later in this paper. 

Very recently there have been several reports of bianisotropic materials which may possess structural 
chirality and exhibit a "negative index" [18, 19, 20, 21], including a commentary on this topic [22]. The 
"negative index" relates to the real part of a wavenumber (relative to that in vacuum), typically correspond- 
ing to propagation in one direction only, for one polarization state only. However, the planewave responses of 
bianisotropic materials is not at all adequately represented that simply: different directions of propagation, 
different polarization states, and the relationship to the time-averaged Poynting vector need to be considered 
too. In contrast, we note that there are some studies in which bianisotropic aspects are taken into account 
better [23, 24, 25]. In the remainder of this paper, we highlight the complications that can arise in the 
planewave response of such a bianisotropic material; furthermore, we demonstrate that some of these com- 
plications can actually be exhibited by relatively simple materials, such as an isotropic dielectric material. In 
so doing, we report on important distinctions between negative refraction, NPV, and counterposition which 
have not been been appreciated hitherto. 

2 Plane waves in a bianisotropic material 

Let us consider the Lorentz-reciprocal [26] bianisotropic material described by the constitutive relations 



with e and /z being the permittivity and permeability of vacuum, respectively, and c = l/^/e /x - These 
particular constitutive relations were chosen because they have been used to describe a material assembled 
from layers of split-ring resonators [27]. This general configuration is a popular one within the negative 
refraction community [21, 20, 25], but it origins predate the current surge of interest in negative refraction 
[28]. The term pseudochiral omega material may be used to describe this material [29]. 

Suppose that a material described by (1) and (2) occupies the half-space z > 0, while the half-space 
z < is a vacuum. We confine ourselves to propagation in the xz plane. In the half-space z < 0, a plane 
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wave with field phasors 



E(r) = E_o ex P [*^o (% sin ip + z cos ip)] 
£L(z) — £Lo exp [ik (x sin + z cos ^)] 



(3) 



is incident on the interface z = 0, where the free-space wavenumber fc = uj^Je a ^ 0l with w being the angular 
frequency. As the incident plane wave transports energy towards the interface, the angle ip G [0, n/2) so that 
the real-valued scalar 

k = k simp e [0, fc„) . (4) 



Two refracted plane waves must exist in the half-space z > 0. Let us represent these plane waves by the 
phasors 

' (7 = 1,2), (5) 



Mr) = Ej exp (ifcj 

= Kj exp (ifc j • r) 



wherein the wavevectors 



(6) 



with € C in general. Thus, the plane waves in the half-space z > arc generally nonuniform. The 
scalars k z j are found by combining the constitutive relations (1) and (2) and the planewave phasors (5) with 
the Maxwell curl postulates. Thereby, we find [10] 



L'E j = 0, 



17 = 1,2), 



with the dyadic 



L = (kj x I + w£ ) • fi 1 • (fc j- x I - wc) + ^ 2 |- 
The dispersion relation det L = yields the two wavenumbers 



(7) 
(8) 



kzl — k \ / l^x I €y 



k Z 2 = k A 



l^zkg 



> . 



(9) 



In these two relations, the square roots must be evaluated such that both refracted plane waves transport 
energy away from the interface z = in the half-space z > 0. 

In order to establish the energy flow associated with the refraction wavevectors kj, the time-averaged 
Poynting vectors 



exp[-2Im (kj) • r] Re . 
• (**x £*-<•£*)]}, 



(i = 1,2), 



(10) 



have to be considered. Both P_\ and P_ 2 lie m the xz plane; furthermore, E_ 1 is directed along the y axis 
whereas E_ 2 lies in the xz plane. 

Let us now consider a specific numerical example in which the constitutive parameters for the bianisotropic 
material occupying the half-space z > are: e x = 0.1 + 0.03i, e y = 0.14 + 0.02z, e z = 0.13 + 0.07z; 
fi x = -0.29 + 0.09i, fly = -0.18 + 0.03z, fi z = -0.17 + 0.6i; and £ = 0.11 + 0.05i. These particular values 
of the constitutive parameters are chosen in order to highlight the complexity of planewave response that 
can be exhibited by bianisotropic materials. Since these materials are artificially constructed materials, 
the range of values that their constitutive parameters can adopt is vast. There is no theoretical barrier to 
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Figure 1: (top) Real part of k Z j/k and (bottom) the angle (in degree) between the real part of k j and the 
positive z axis, both plotted against n/k Q . The solid curves correspond to j = 1, and the dashed curves to 
j = 2. As Ymkzj > for all n/k £ [0, 1) and both values of j, the imaginary parts of k z i/k and k z ijk a 
have not been plotted here. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 



0.0 0.1 0.2 0.3 0.4 

x/k 



0.5 0.6 



Figure 2: (top) The x component of the normalized time-averaged Poynting vector and (bottom) the angle 
(in degree) between the time-averaged Poynting vector and the positive z axis, both plotted against n/k . 
The solid curves correspond to j = 1, and the dashed curves to j = 2. As | • Pj > for all n/k £ [0, 1) 
and both values of j, those quantities have not been plotted here. 



the particular values chosen here: these describe a dissipative pseudochiral omega material. It is through 
accessing unconventional values of the constitutive parameters that we determine if mctamaterials may 
exhibit their exotic, and potentially useful, properties. 

In Fig. 1, the real part of k Z j/k , (j — 1, 2), is plotted as a function of (/t/fe ) £ [0, 1). As the imaginary 
parts of both k z \ and k Z 2 turned out to be positive, both refracted plane waves must attenuate as z — > oo, 
in consonance with our understanding of a passive medium. Also plotted in Fig. 1 is the angle between the 
real part of kj, (J — 1, 2), and the positive z axis. The refracted plane wave labeled 1 is positively refracted 
for < (n/kn) < 0.14 but negatively refracted for 0.14 < (n/ko) < 1. Additionally, the refracted plane wave 
labeled 2 is negatively refracted for < (n/k ) < 0.22 but positively refracted for 0.22 < (n/k ) < 1. 

The normalized x component of the time-averaged Poynting vectors for both refracted plane waves are 
plotted against {n/k ) in Fig. 2. The angle between (j = 1, 2), and the positive z axis is also plotted. The 
z components of Pj and P_ 2 are positive for all (n/k ) £ [0, 1), in accordance with the rule to evaluate the 
square roots in (9). But, whereas the x component of P 1 is negative for all (n/k ) £ [0, 1), the x component 
of P 2 is negative only for (n/ko) £ [0,0.01). Therefore, P_ x always subtends a negative angle to the positive 
z axis whereas the sign of the angle that P 2 subtends depends on k. 

The phenomenon of counterposition [16] arises when the inequality 

Re (k zj ) x • Pj < (11) 

is satisfied. By comparing Figs. 1 and 2, we see that counterposition occurs for the refracted plane wave 
labeled 1 when < {n/k ) < 0.14. That is, the K-range for counterposition coincides with the K-range for 
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Figure 3: The quantity Re (fc^) • P- / (fc | |) plotted against {n/k ) £ [0, 1). The solid curve corresponds 
to j = 1, and the dashed curve to j — 2. 

positive refraction. However, this is not the case for the refracted plane wave labeled 2: here, counterposition 
occurs only for 0.01 < (re/fc„) < 0.22. 

The quantity Re {k j) • P - determines whether the phase velocity of the refracted plane wave labeled 
j is positive or negative. In Fig. 3, the quantities Re (fc ) • P- / (fc n | P^ |), (j = 1, 2), are plotted against 
(n/k ) £ [0, 1). The phase velocity is negative for the refracted plane wave labeled 1 for 0.07 < {n/k a ) < 1, 
while for the refracted plane wave labeled 2 it is negative for < (n/k Q ) < 0.04. Thus, the K-ranges for NPV 
do not coincide exactly with those for negative refraction. For example, consider k = 0.1fc o : (i) The angle 
between Rek 1 and the +z axis is 86° while the angle between P 1 and the +z axis is —13°. This refracted 
plane wave has NPV but it is positively refracted, (ii) The angle between Refc 2 an d the +z axis is —79° 
while the angle between P 2 an d the +z axis is 58°. This plane wave has positive phase velocity but it is 
negatively refracted. 

Our results are summarized in Tables 1 and 2 for the plane waves labeled 1 and 2, respectively. These 
tables include illustrations of the directions of Rekj and P •, (j = 1,2), for representative values of k. In 
the Appendix we have confirmed that the corresponding reflection and transmission coefficients are nonzero. 

3 A simple case: an isotropic dielectric material 

An important point demonstrated in the preceding section is that, under certain circumstances, NPV and 
positive refraction can co-exist for the bianisotropic material under consideration. This co-existence of NPV 
and positive refraction has previously been predicted for certain photonic crystals [30] and nondissipative 
uniaxial dielectric mediums described by indefinite permittivity tensors [14] (as well as for certain active 
materials [31]). In fact, the independence of NPV and negative refraction is also manifest in simpler materials. 
Let us take, for example, the simple case where the half-space z > is occupied by an isotropic, dielectric 
material with relative permittivity scalar e. This represents the simplest specialization of (1) and (2). Here 
there is only one refraction wave vector, i.e., k 1 = fc 2 , where \k 12 \ = k y/e, but the two associated time- 
averaged Poynting vectors, namely P 1 and P 2 , are generally distinct in terms of both magnitude and 
direction for e £ C. As previously, P_ 1 corresponds to E 1 being perpendicular to the plane of incidence and 
P 2 corresponds to E_ 2 being parallel to the plane of incidence. 

For the purposes of illustration, let us take the relative permittivity scalar e = — 6+2. 5i. A straightforward 
calculation reveals that the refraction is always positive for both perpendicular and parallel polarization 
states. Furthermore, when the electric field phasor is perpendicular to the plane of incidence, the phase 
velocity is positive and there is no counterposition, regardless of the value of Indeed, the time-averaged 
Poyting vector and the phase velocity are parallel for this polarization state. A different picture emerges 
when the electric field phasor is parallel to the plane of incidence: counterposition occurs for all ip and 
the phase velocity changes from positive to negative as tp increases in value, with the phase velocity being 
orthogonal to the time-averaged Poynting vector when i/j = 32°. 

The directions of the refracted wavevectors and the corresponding time-averaged Poynting vectors are 
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Table 1: Ranges of values of k for which the plane wave labeled 1 supports negative/positive refraction, 
counterposition and negative/positive/orthogonal phase velocity. The final column shows the directions in 
the z > half-space of Rek 1 (thick dashed arrows; blue in electronic version) P_ l (thick solid arrows; red 
in electronic version) for representative values of k; also shown are the directions of the wavevectors for the 
incident and reflected plane waves in the z < half-space (thin solid arrows; green in electronic version). 
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Table 2: As Table 1 but for the plane wave labeled 2. 
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Table 3: As Table 1 but the z > half-space is replaced by an isotropic dielectric material with a relative 
permittivity e = —6 + 2.5i. The plane wave is polarized perpendicular to the plane of incidence. Notice that 
here Re fc 1 and P_\ are parallel. 
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Table 4: As Table 3 but the plane wave is polarized parallel to the plane of incidence, 
illustrated in Tables 3 and 4 for both polarization states and for some representative values of k. 

4 Concluding remarks 

The combination of anisotropy and magnetoelectric coupling can result in a much more complicated planewave 
response than is associated with isotropic dielectric-magnetic materials. This is especially significant when 
exotic constitutive parameters ranges — such as those associated with certain materials which support neg- 
ative refraction — are considered. A further level of complication is introduced by considering nonuniform 
plane waves. The relationships among the three phenomcnons of negative refraction, NPV, and counterpo- 
sition highlight the complications. In particular, our investigation based on a pseudochiral omega material 
has revealed that: (a) negative refraction can arise even though the phase velocity is positive, and positive 
refraction can arise even though the phase velocity is negative; (b) counterposition can arise in instances of 
positive and negative refraction; (c) whether or not positive or negative refraction arises can depend upon 
the angle of incidence; and (d) at the transition from positive to negative phase velocity with increasing 
angle of incidence, the phase velocity and time-averaged Poynting vector are orthogonal to each other. 

The exhibition of NPV in an isotropic dielectric material — which follows due to the consideration of 
nonuniform plane waves — is particularly noteworthy. If only uniform plane waves are considered then this 
outcome is impossible: magnetic properties or anisotropy, for example, would also be needed in order to 
support NPV [12, 14]. 

The findings reported herein further emphasize the importance of fully characterizing bianisotropic ma- 
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tcrials, instead of attempting to do so with a single scalar refractive index. Also, the complications that 
can be introduced by taking account of nonuniform plane waves, even in isotropic dielectric materials, are 
highlighted. The demonstration that NPV and negative refraction can arise independently of each other 
has highly significant consequences for researchers exploring the realm of materials which support negative 
refraction, bianisotropic or otherwise, and beyond. 
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Appendix 



We present here the reflection and transmission coefficients corresponding to the reflection-transmission 
problem. 

In terms of linearly polarized states, the incident plane wave is represented by the electric and magnetic 
phasors 

ELiiz) = (ass + o-pp^j cxp [i (kx + k a z cos ip)] 

1 f ~ v \, z<0, (12) 

H_i(l) — — \ a sP + ~ a pS-j ex P ( KX + k z cos ip)] 

where ry = ^/^ /e and the unit vectors 



s = y, p — =F xcosip + z_smip. 
The corresponding reflected plane wave is represented by 

E r (r) = (r s s + r p p_^j exp [i (kx — k zcosi/))] 

H_ r (r_) = — ( r sP_ — r p s) cxp [i (kx - k z cos ip)] 

and the transmitted fields are written as 

E-tiL) — t\s exp [i (kx + k z \z)\ + t 2 p exp [i (kx + k Z 2z)] 

S_ t (r) = txfjT 1 • ^— ki x L-(^J ' 1 exp [i (kx + k zl z)] 



+^2M 1 • ( — k 2 x J — ( ) • p cxp [i (kx + k Z 2Z 



wherein vector 
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satisfies p • p* = 1. 

The amplitudes r s p and tip of the reflected and transmitted plane waves are conveniently related to the 
amplitudes a s , p of the incident plane wave via the matrix relations 

(17) 
(18) 

In order to find the reflection and transmission coefficients (r ss etc. and ti s etc.), the boundary conditions 
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are imposed. Thus, the four simultaneous linear algebraic equations 
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(20) 
(21) 

i42Q) 
(23) 



emerge, which may be solved to provide r ss etc. and ti s etc. by straightforward algebraic manipulations. 



For the bianisotropic metamaterial considered here, it transpires that r, 



sp 



' ps 



and tip = t2. s 



The absolute values of : 



0. 

ss) r ppi tis and t2p are plotted against n/k in Fig 4. The transmission coefficients 
have magnitudes > 1 in the chosen representation of the transmitted fields, but we have verified that the 
principle of conservation of energy is satisfied by the results. 
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